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Abstract 

We study nonlinear cosmological perturbations during the post-inflationary evolution, using the 
equivalence between a perfect barotropic fluid and a derivatively coupled scalar field with Lagrangian 
[— (d4>) 2 ]( 1+w ^ 2w . Since this Lagrangian is just a special case of ^-inflation, this approach is analogous 
to the one employed in the study of non-Gaussianities from inflation. We use this method to derive 
the second order metric during matter dominance in the comoving gauge directly as function of the 
primordial inflationary perturbation £. Going to Poisson gauge, we recover the metric previously 
derived in the literature. 



1 Introduction 

Recently, there has been an extraordinary improvement in the accuracy of cosmological observations, 
especially regarding the statistical properties and the evolution of perturbations around a FLRW 
Universe. In order to fully exploit the data, it has become necessary to go beyond linear perturbation 
theory. This is mandatory, for instance, if one wants to study primordial non-Gaussianities which, 
by definition, are sensitive to the non-linear interactions of the theory. 

The calculation of primordial non-Gaussianities in various models is based on a perturbative 
expansion of the action - or of the equations of motion - of one or more scalar fields. In general, the 
mixing of these scalar fields with gravity cannot be neglected. This can be done by using an ADM 
approach which allows to write an action for the relevant degrees of freedom taking into account 
gravitational perturbations pp. In this paper we want to extend this formalism to the study of 
cosmological perturbations after inflation, during the standard FLRW evolution. 

During its evolution, the Universe is filled with cosmological fluids such as radiation or cold dark 
matter. At first sight a fluid is very different from a scalar field, so that it is not clear how to extend 
the aforementioned approach to our case. However, in most cases the cosmological fluids can be 
considered as perfect and barotropic. In this case, as we will review in section [21 the dynamics of the 
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fluid is exactly equivalent to the one of a derivatively coupled scalar field, playing the role of the fluid 
potential velocity. Using this equivalence, the usual treatment of perturbations in the primordial 
Universe, used when they are generated and stretched out of the Hubble radius, can be extended to 
their subsequent evolution in the FLRW Universe, when the modes re-enter the Hubble radius. 

In this article this approach is applied to find the exact metric at second order during the matter 
dominated era. Our results are compatible with those obtained previously in [3], by directly 
working at the level of the Einstein equations. The knowledge of the second order metric is the first 
step for the calculation of any observable beyond the linear approximation. This calculations are 
quite relevant: in the case of the CMB anisotropy the expected magnitude of a generic second order 
effect is comparable to the sensitivity of the forthcoming Planck satellite. 

As we will see, the calculation of the second order metric will closely parallel the one of the 
three-point function of scalar and tensor perturbations generated during inflation. In particular, we 
will make use of the action calculated for this purpose for fc-inflation [H [5], expanded up to third 
order. In this context, it is natural to use the comoving curvature perturbation £ as the variable 
describing the scalar perturbations, the same which is commonly employed for the quantization of 
scalar fluctuations during inflation. This variable is also well-known to be nonlinearly conserved on 
super-Hubble scales independently of the details of the cosmological evolution [6]. For this reason it 
is the natural variable to set the initial conditions for the post-inflationary dynamics. Therefore, an 
advantage of our approach is that the dynamics of perturbations is described by the same variable 
and the same approach, from its generation at early times to the late times observations. 

As a warm-up, we start with linear perturbation theory in section[3l Then, in section^ we extend 
our calculation to second order for scalar and tensor modes, assuming that primordial gravitational 
waves are negligible. In this way we obtain the full second order metric in the comoving gauge 
including tensor modes generated by scalars. In section[5]we check our results with the ones obtained 
in Poisson gauge (2j [3] after a suitable second order gauge transformation. Conclusions are drawn 
in section [6j 



2 A fluid as a scalar field 

In cosmology, when dissipative phenomena are negligible, the energy content of the Universe can be 
approximated as a sum of perfect fluids. A perfect fluid is defined to have a stress-energy tensor of 
the form 

Tfw = (p + p)u^u u + pg^v , (1) 

where p and p are the energy density and the pressure, while is the fluid 4- velocity. Cosmological 
fluids are usually taken to be irrotational. This assumption is justified by the absence of vorticity 
in the initial conditions set by inflation and by the fact that vorticity is diluted by the expansion 
of the Universe. The fluids are also taken to be barotropic, i.e., their pressure is a function of the 
energy density only, p = p(p). Under these conditions, each fluid is characterized by a single scalar 
function, so that it is not surprising that its dynamics can be described in terms of a scalar field. 

Indeed, let us consider a derivatively coupled scalar eft in Minkowski spacetime with Lagrangian 
densit}0 

c = p(x) , x = -a^m . (2) 



1 For Lagrangian approaches which describe also vorticous motion and non-barotropic perfect fluids see 
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Varying the action yields the equation of motion 



d fi [P'{X)d^] = . (3) 
The stress-energy tensor of this field is given by 

= 2P'(X)d^d u ct> + P(X)g^ , (4) 
which is of the perfect fluid form ([1]) if we identify 

p = 2P'X-P, p = P, = (5) 

VA 



The perfect fluid interpretation makes sense only if d^cf) is everywhere timelike and future directed. 
Projecting the conservation equation of the stress-energy tensor of a fluid, d^Ty = 0, along and 
orthogonal to the fluid flux yields the equation of conservation of energy and the Euler equation. In 
the case of the stress-energy tensor ^ , the Euler equation is a trivial identity, while the conservation 
of the energy is equivalent to the equation of motion ([3]) . 

The equation of motion ([3]) can be interpreted as the conservation of the current 

J» = 2^XP'(X) ■ . (6) 

This conservation is a consequence (by Nother theorem) of the invariance of the action under shift of 
<j> '. <j> —* + const. From the fluid point of view this current describes the conserved particle density 
flux = mi' 1 , where n is the number particle density. Therefore one can identify n = 2yXP'(X). 
This yields a physical interpretation of the norm of d^cp: yX = (p + P)/n is the so called specific 
inertial mass [7, 8j. It is also straightforward to verify that 

so that yX acts as a sort of conjugate variable with respect to n [9j. 

It is well known that for a perfect fluid the entropy per particle is conserved along the fluid flow. 
This can be checked using the continuity equation d^J^ = 0, i.e. eq. ((5J, and the energy conservation 
u^d v T^ = [12] . Furthermore, as we have discussed, our approach describes a barotropic fluid which 
implies that the entropy per particle is everywhere constant. In other words, the Lagrangian (|2|) can 
only describe mechanical excitations of the fluid. It cannot take into account dissipative irreversible 
processes like heat conduction or viscosity. 

We will be interested in studying perturbations around a homogeneous configuration. In Min- 
kowski spacetime (j) = ct is a solution of the equation of motion for any c. Different values of c 
describe different unperturbed values of the energy density p(c 2 ). The dynamics of fluctuations 
around this background can be studied by expanding the Lagrangian using eft = ct + 6cf)(t, x), where 
8<p describes the compressional mode of the fluid. At second order we obtain 



C = P'(c 2 )[6<f> 2 - (V5cj)) 2 } + 2P"(c 2 )c 2 5<J) 2 . (8) 



13 El El- 
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From this expression one sees that the speed of sound of the excitations is given by 



P'(X) + 2XP"(X) 

which, as expected, is the usual adiabatic speed of sound in a barotropic fluid 

2 V'(X) 



(9) 

X=c 2 



= t ■ do) 

X=C 2 d P 



Note that although the Lagrangian ([2]) is Lorentz invariant, this symmetry is spontaneously broken 
by the vacuum <j> = ct. For this reason the speed of perturbations will in general differ from the 
speed of light cr s = 1. 

A standard case of barotropic fluid is given by a linear equation of state p = wp with w = const. 
From the first two equalities in eq. (J5j) one deduces (see for instance [TUl El] ) 

P = X 1 ^ , w / , (11) 

up to a proportionality constant which is irrelevant for the classical theory. In this case, the speed of 
sound in eq. (|10p reduces to c 2 = w. As an example, one can consider a radiation fluid with equation 
of state w = 1/3. In this case the Lagrangian (jlip reduces to 

c = x 2 = {-d^d^f . (12) 

The inclusion of gravity is completely straightforward. The equation of motion ([3]) becomes 

dpW=jP , {X)&><j>) = . (13) 
In an expanding FLRW background, the homogeneous solution satisfies now 

n = P'4> oc a -3 . (14) 
Indeed, for a constant w this is the standard redshift of the energy density, 

p ex ^ ex a ~ 3 ^ . (15) 



To summarize, we are able to describe the dynamics of one or more fluids coupled with gravity, 
within a Lagrangian formalism. This can be employed to study cosmological perturbations in the 
presence of barotropic fluids. Note that the Lagrangian ([2]) is a particular case of the so called 
k— essence/A;— inflation scenarios |13} I14j. For this class of Lagrangians, cosmological perturbations 
have been extensively studied starting from [TT] and more recently extended to second [HE] and third 
order [15J to study the non-Gaussianities produced during inflation. In this paper these calculations 
can be reinterpreted as describing the nonlinearities of a fluid. 

In the following we will be interested in the study of dark matter perturbations, which would 
correspond to the dust case w = 0. Obviously the zero pressure limit must be taken with care, as 
from eq. © the Lagrangian strictly vanishes for a pressureless fluid. 
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3 Background and linear dynamics 

Following the discussion of the last section, the dynamics of a perfect fluid coupled with gravity can 
be described by the action 

S = ^J a 4 x^-g[R + 2P{X)\ , (16) 

where we chose units such that Mp 2 = 8itGn = 1. In particular we are interested in the case p = wp 
with constant to, when P is given by eq. (jllj) . 

We assume a flat FLRW metric ds 2 = —At 2 + a 2 (t)dx 2 . Friedmann equations read 

H 2 = - (2XP' - P) , (17) 
3 

2H + m 2 = -P , (18) 

where H = a/a is the Hubble expansion parameter. It is useful to define 

H 3 . „ . . 

e ^-Jp = 2 ( - 1 + w) ■ (19) 
Note that, although the notation is inspired by the inflationary case, we are not assuming that e is 
small H 

Following Maldacena pQ, the study of perturbations can be done using the ADM splitting of the 
metric [IB] . 

ds 2 = -N 2 dt 2 + h ij (dx i + N i dt)(dx j + N j dt) . (22) 

The action (|16p becomes 

S = i J dtd 3 xVh N(R( 3) + 2P) + iV" 1 - E 2 ) , (23) 

where R^> is the curvature scalar computed from hij, E = E\ and 

/••';, ^f/'; ; V..V, Yj.X.) . (24) 



2 

The covariant derivatives Vj are with respect to the 3-metric hij and all roman indeces . . . are 
raised and lowered with this metric. 

The action (|16p describes 3 dynamical degrees of freedom: one scalar mode for the fluid exci- 
tations and the 2 tensor helicities of the gravity waves. In the ADM formalism, these degrees of 
freedom are contained in the scalar field <j) and in the 3-metric hij, while the lapse N and the shift 
Ni are not dynamical. Following [I] we choose the gauge 

5(j) = 0, hij = a 2 e 2C -hij , hij = Sij + 7^ + §7i/7/j + ■■■ , 

(25) 

det h = 1 , ju = , dtfij = , 



2 To simplify the comparison with the literature we give here the parameters used in refs. [11 [5] for the w = 
const case: 

,= -L=0, «=1 1=1-1, s =If£ =0 , (20) 
eH c% w H c s 

and 

E = XP' + 2X 2 P" = — , \ = X 2 P" + -X 3 P"' = -— (l- -) . (21) 
w 3 6w V w I 
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where £ describes the scalar mode and 7 the tensor ones. This gauge is comoving in the sense that 
the fluid 4-velocity is everywhere orthogonal to the equal time surface^. 

As the action does not contain time derivatives of N and N% , these variables act as Lagrange mul- 
tipliers, i.e. their equations of motion are algebraic constraints. These equations are the momentum 
and Hamiltonian constraints of Einstein equations: 

V, [N" 1 (E) - 8}E)] = , (26) 

#(3) + 2P - 4XP' - (E ij E ij - E 2 ) = . (27) 

To warm up let us calculate the action for the scalar mode Q at second order, in order to find 
the linear equation of motion and the first order metric. As there is no mixing between Q and 7 at 
second order in the action, in this section we can set 7 = 0. To compute the second order action, we 
need to solve the constraints (|26l27p and plug their solution for iV and iVj back into the action (|23p . 
This needs to be done at first order only, since the second-order solutions for N and N l will multiply 
8C/5N or 5C/5N 1 at zeroth order, which vanish on the background [lj. To solve the momentum 
constraint (|26p we decompose iV, as N{ = diip + Nt% where diNxi = 0. By defining N = 1 + 5N, one 
finds, at first order, 

8N = jj, N Ti = 0. (28) 
Furthermore, one can find ip at first order by solving the energy constraint^ (j27j) . 

i> = -! + — d~H ■ (29) 
H w 

Note that the last term of the previous equation contains w at the denominator so that one is not 
allowed to take the limit w — > at this stage. 

Substituting the solutions (|28l29p into the action (|23|) one obtains, after some integration by 
parts, the second order action for £, 



S'o - / dtd 3 xa 3 - 
w 



c--(dcr ■ (30) 



The equation of motion derived from this action is thus 
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in 

( + 3H(--d 2 ( = 0. (31) 



We are interested in the limit w — > 0. Setting w = 0, the growing solution of eq. (13ip is simply 
a constant, £ = £0, where Co is the perturbation generated during inflation, which remains constant 
on super-Hubble scales independently of the equation of stat^E However, to find the metric we need 



3 We follow the notation of pQ using £ to denote the curvature perturbation in the comoving gauge. Often 
this quantity is instead called 1Z reserving £ to denote the curvature perturbation in the uniform density 
gauge. 

4 Here and in the following d 2 = didi and d~ 2 is its inverse. 
5 For a generic equation of state p(p) the action for £ is 



S = I <\td 3 x a 3 4 



c 2 - i(dc? 

a 1 



(32) 
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to plug the solution for ( into eqs. (|28l29p and for £ in eq. ([29]) we need to keep the leading order 
correction in w in eq. (|3ip . 

The general solution of this equation in Fourier space can be written as a linear combination of 
Bessel functions, 

Cg(r) = dr-^MkrV^) + C 2 r^Y v {kr^), v = ^+sw) ' (34) 

where we have introduced the conformal time r = J dt/a. The growing mode of this solution goes 
as 

( = (o + in ^d 2 ( + 0(w 2 ), (35) 
where Co is the asymptotic value of ( on super-Hubble scales. The time derivative of (f35l) is 

C = ^ d2 Co + 0{w 2 ) . (36) 
This can be plugged into eqs. (|28l29p which yield, at lowest order in w — > 0, 

SN = 0, 1> = ~\%- (37) 
5 H 

The metric at first order in perturbations is thus given by 

ds 2 = -dt 2 - -^-diCo dtdx i + a 2 (l + 2( )dx 2 . (38) 
5H 

Note that under a coordinate transformation 

t^t+TS' x^x*, (39) 
bid 

this metric takes the Newtonian gauge form 

ds 2 = -(l + 2$)dt 2 +a 2 (l-2^)dx 2 , (40) 

with 

$ = * = --Co, (41) 
5 

which is the well known linear relation between the Newtonian potential and Co iu matter dominance. 



where c 2 s is given by eq. ([9]). This gives the following equation of motion 

f+<sH*g)-;J* , <-°- < 33 > 

From this equation we see that £ is constant on large scales, k/aH <C 1, independently of the energy content. 
This is not true, for instance, for the Newtonian potential. 
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4 Second order perturbations 



4.1 Scalar perturbations 

We are interested in deriving the third order action for £. In order to do that we follow the same 
procedure we did at linear order and initially set 7 = everywhere. Tensor modes generated at 
second order by scalar perturbations will be considered in the next section. To derive the third order 
action we only need the solutions to the constraint equations at linear order, i.e. eqs. (|28|29p . This 
is because the second and third order solutions would multiply the first and zeroth order constraint 
equations, respectively p]. We will be interested in the second order solutions of the constraint 
equations only later, when we will derive the explicit expression of the metric at second order. 
It is straightforward to derive 



E ij E ij - E 2 



\2 - e 



6(H + cr + 4^r(H + C) (diNi + di(Ni 



-K 



-(diNj + djNiY - (diNj + djNiXNidjC + NjdiQ - (diNiY 



and 



i?( 3 ) 



-2^[2d 2 ( + (d() 2 ]. 



(42) 
(43) 



Expanding P(X) with X = cfr/N and using the first equality in (|28p one obtains [I] 



P(X) 



P + P' 



c c 2 c 3 c 



X 



+ \P" (4^+8^4-12^ 



X< 



(44) 



where, on the right hand side, X, P and its derivatives are evaluated at zeroth order. The lapse has 
been split as N = 1 + 5Ni + 5N2, where 5N\ = (/H is the first order contribution, eq. (|28p . while 
<5iV2 is the second order one. As we discussed, terms containing 5N2 will not appear in the action; 
however they will be important to derive the metric. 

Plugging the previous expressions into the action and using the solutions to the constraint 
equation (I28p one obtains, after several integrations by parts, the third order action for 



dtd^xa- 5 



w 



(2 + i)C 3 -4(<90 2 C 



3Hw 



w 



—diCd^d 2 ^ + {d t d^d t d^ - (d 2 ip) 2 ) (3C - 4) 



H' 



(45) 



(46) 



where ip is given by eq. 

^ = -77 + — d- 2 (. 
H w 

This action describes the cosmological scalar perturbations of a perfect fluid with constant equation 
of state w. Notice that it is just a particular case of the general action derived in refs. [U [5] in 
the context of generalized inflationary models, although it describes here a quite different physical 
system. 



S 



The action is simplified by a field redefinition, Cn = Q — /(C), where [5] 



+ 



' aiCaidr^-a^aidjiaiCdjd- 2 ^ . (47) 



2# 2 



As / is quadratic in C, the field redefinition does not change the second order action. Note that the 
function /(C) contains either spatial gradients or time derivatives of C- Thus, for any constant w, 
the field redefinition (|47p vanishes on super-Hubble scales, so that in this regime C and Cn coincide. 
In terms of Cn the action reads 



/ 



dt d 3 xa 3 - 
w 



9 ( - " 1 ) I + I ( 3 " - ] Ul + t^(5 + u;)C„(flkC "' 2 



w 



2 - - - Cn^C^^Cn + -rd'Udtd-'Cn 

' ' w Aw 



II ' 2 V" «•>/ s " x " 2a2 v " ' 

(48) 



We are interested in a dust fluid w = 0. The equation of motion at second order derived from 
this action can be solved perturbatively by plugging the first order solutions C = Cn = Co into the 
second order terms. At lowest order in w — * one finds 

Cn + 3#Cn = , (49) 

where 

F = -Co<9 2 Co - ttt(^Co) 2 - ^cT^Co^Co) • (50) 
lb » 



The right hand side of eq. (|49p is negligible at very early times, when all the modes are well out of 
the Hubble radius, k/aH <C 1, so that Cn = const is a solution. As in this regime C and Cn coincide as 
discussed, this corresponds to the well-known fact that C is conserved on super-Hubble scales, even 
nonlinearly. In this regime the Fourier modes of C are decoupled and equal to the initial condition 
set by inflation, C = Co for k/aH <C 1. Solving the equation above and going back to £ using (|4T|) 
yields 

C = Co - ^jpd^d^iW.Co) ■ (51) 

This expression gives C at second order during matter dominance, as a function of its initial condition 
Co- We postpone its discussion to the end of this section and we proceed to conclude the calculation 
of the second order metric: we need to solve for and Ni at second order. 

As at first order, to solve for the whole metric we need to keep terms of order w in the second 
order evolution equation which now reads 

( n + 3H( n -^d 2 Cn = ^(F + wF 1 ) + ^G , (52) 



where 



F x = -^d 2 C Q + ^{d l ^f + \d i d- 2 {d 2 C d l Co) , (53) 

G = -§(9 2 Co) 2 -|Co5 4 Co + ^W^ 2 Co 
25 5 20 

+ m d2{diCo)2 ~ 4^" 2(w4Co) • (54) 



9 



To obtain the w corrections to the right-hand side we used the linear evolution of £ at order w 2 : 

2w 



C 



(55) 



Solving eq. 

c 



5a 2 H su 25a 2 H su 35a 4 # 3 
and writing the solution in terms of £ using the field redefinition (|47|) one finds 



Co 



+ 



1 



d 2 didjidiCodj Co) 



5a 2 H 2 

2d 2 (o - 4Co5 2 Co + ^Co) 2 + ^^0- 2 (a 2 Cockn ! 

5 



5a 2 F 2 



+-^- 2 (5,C 5 4 Co) + ^diididjCodj&Co) 



(56) 



Now that we have a solution for £ including O(w) corrections, we can proceed to solve the 
constraints to get the final second order metric. The momentum constraint is a vector equation; 
taking its transverse part with the projector 6{j — didj/d 2 we obtain an equation for the transverse 
part of the shift vector Nxi {diNxi = 0), 



N Ti = ~^d- 2 [W 2 ^ ~ d- 2 did k (d k ( d 2 C )] 



(57) 



Notice that the transverse part of the shift vanished at first order, see eq. (J28J). The longitudinal 
part of the momentum constraint gives an equation for the lapse function at second order, 

^ [(d(o) 2 ~ 4Co5 2 Co] 



5N 2 



5a 2 



2w 



[3(d 2 Co) 2 + mCodid 2 ^ + mdjCof] + 0{ 



w 



(58) 



175a 4 tf 4 

Although the lapse perturbation vanishes in the limit uu — > 0, its expression at order w is necessary 
to solve the energy constraint, similarly to what happened at first order. This gives the second order 
correction to the shift function 

1 



1p 2 



5H 



-2 



(di( y -sd-'didjidiCodjCo) 



+ ^4^(l(d 2 Co) 2 + W,d 2 Co + y(a^Co) 2 



(59) 



5a 2 # 2 V7 

Plugging in the ADM metric (|22p the first order results (|28p . (|29p and the second order ones 
56|) . ([57|) . (|58|) and ([59]) . we finally obtain the second order metric for w = 0, 

4 



500 

9oi 



-1 + 
1 

~5H' 



25a 2 H 2 



(5,Co) 2 



(60) 



2Co - d- 2 (djCo) 2 + Sd^djdkidjCodkCo) 

-^d- 2 (l(d 2 c ) 2 + W t d 2 Co + ^ 



4 1 



-^d- 2 [W 2 ^ - d- 2 dMdkCod 2 Co)] , 



9i.i 



1 + 2Co + 2C 2 



5a 2 F 2 



d 2 d k di(dk(odiCo) 



a 2 lij 



(61) 
(62) 
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Before moving to the calculation of the tensor contribution 7, let us discuss eq. (|5ip which 
describes the second order evolution of £. In Fourier space this equation reads 




d 3 q (k ■ q)(k 2 — k ■ q) 




(63) 



* 5a 2 H 2 J (2vr) 3 




It is important to stress that the kernel inside the integral does not vanish in the limit k — > 0. 
Surprisingly, this implies that a very long wavelength mode well outside the Hubble radius evolves 
in the presence of short wavelength perturbations: £ is not conserved on super-Hubble scales as 
its second order contribution in the matter dominated era grows like the scale factor a oc l/a 2 H 2 . 
However, this is not in contrast with the literature on the conservation of C, [6], as £ is conserved 
when all the modes are out of the Hubble radius, so that the second term on the right hand side of 
eqs. (|5T1) and ([63]) can be neglected. 

The fact that Q is not conserved in our case is not so surprising after all. Indeed, at nonlinear 
order one can define an infinite set of variables which coincide with and which are thus conserved, 
when all the gradients can be neglected, i.e. out of the Hubble radius. However, when sub-Hubble 
modes are considered, if one of these variables remains constant, the others will not in general. In 
fact, one can trivially define a variable which is conserved on all scales, i.e., 



even though this has not a particular geometrical meaning and we do not expect it to remain 
conserved when a different equation of state is considered. In order to discuss the physical significance 
of the variation of Q one needs to compute at second order the relationship with a measurable 
quantity, like Cosmic Microwave Background temperature fluctuations. This goes beyond the scope 
of this article. 

4.2 Gravitational waves 

So far we did not discuss gravitational waves. At linear order tensor and scalar modes are decoupled, 
so that gravitational waves can be completely neglected in the limit where the primordial contribution 
generated by inflation is very small. Beyond the linear approximation scalar and tensor modes mix, 
so that it is not consistent to set 7 = 0. In the following we are going to assume that the amplitude 
of primordial tensor modes is very small: gravitational waves will only be generated by the scalar 
modes through couplings of the form 7CC- 

To study the generation of gravitational waves we need the quadratic action for 7 and the cubic 
terms of the form 7CC The constraint equations will not enter neither in the derivation of the action 
nor in obtaining the expression of the second order metric. Indeed 7 does not appear at first order 
in the constraint equations and quadratic terms are negligible as 7 = 0(C 2 )- 

Let us start with the quadratic action for 7. Making use of the expressions 



1 



d^didjidiCdjQ , 



(64) 



c = c + 



5a 2 H 2 



EijEV -E 2 D 



(65) 



and 




(66) 
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after integration by parts this yields 

Sjj = ~ J dtd 3 xa 3 \ - -^dkjijdkjij] . 



(67) 



For the CCl terms we use the expressions 



R {3) D -^-ijihOjC - ^lijCdidjC 



and 



e ^ N -l a 3 - E 2 ) d -ajij 3( 



H 



didjif) + -dkjijdidjipdkip , 
a 



and we obtain 



S. 



7CC 



dt d 3 x a 3 



Ha 2 



1 



C 



1 



H 



2a 4 



(70) 



This part of the action is the same as the one derived by Maldacena pQ for a scalar field with 
standard kinetic term. This is not surprising as these couplings do not depend on P{X) in (|23|) . 

The equation of motion for 7 can be obtained by varying the actions (|67|) and (|7U|) with respect 
to 7y. Here there are no subtleties in the limit w — > so that one can drop terms containing £. 
Notice that 7^ is constrained to satisfy the transverse and traceless conditions d^ij = and 7^ = 0. 
As such, the 7's in the action above must be understood as projected by the transverse and traceless 
projector, i.e., 

(71) 



P; 



TT 



ij kl 



2 {PikPjl + PjkPil — PijPkl) , 



where Pij is a symmetric transverse projector given by 



Pi, 



Thus, the evolution equation reads 



n- 



TT 

ij kl 



jd k ( di( + 



2ba A H 2 



d 2 (dkCadtto) 



(72) 



(73) 



where we have simplified the term in the square braket on the right hand side, taking into account 
that we are interested only in its transverse component. The solution of this equation, with the 
appropriate initial condition 7^ = 0, reads 



lij 



9 



1 ji(fcr) 
3 kr 



d- 2 + 



1 



ha?H 2 



P?l (dMCo) 



(74) 



where the spherical Bessel function ji(x) can be written in terms of trigonometric functions as 
ji(x) = sin(x)/x 2 — cos(x)/x. 
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5 Transforming to Poisson gauge 

In this section we want to write our second order metric in conformal Poisson gauge, i.e. in the form 
ds 2 = a 2 (r) {-(1 + 2$ P )dr 2 + 2u^ i drdx i + [(1 - 2* P )(5y + j Pij ] dxMx J '} , (75) 

where ujp { is transverse, diiOp i = 0, and 7p is transverse and traceless. Going to this gauge will 
enable us to compare our results with those obtained in [2j [3] . 

For the second order gauge transformation we use the notation of [2] (apart from exchanging the 
name of the scalar potentials $ and and we write our metric as 

ds 2 = a 2 (r) {-(1 + 2$ ? )dr 2 + 2(a/^ + d i u < )dTdx i + [(1 - 2* c )<^ + 7y ] dxMx j } , (76) 

where 



r(^Co) 2 
1 



25a 2 H 2 

2 Co + T^d- 2 [(diCo) 2 - sd^didjidiCodjCo) 



5aH 
4 



25a 3 # 3 
4 



5aH 

d- 2 



^(d 2 c ) 2 + d t ( d t d 2 Co + ^(^Co) 2 



5a# 



a- 2 [^Co5 2 Co - d-tdidjidjto&to)] , 



-Co - Co 2 + 



i 



5a 2 H 2 



d^didjidiCodjCo) , 



(77) 

(78) 

(79) 
(80) 



and jij is given by eq. (|74p . The gauge transformation up to second order is, in conformal coordinates 



(r,x l ), 



(81) 



where at each order £° = a and = d l (3 + cf with the vector d l transverse, didi = 0. The first order 
piece is fixed by eq. ([39]) . i.e. a\ = —2(o/(5aH) and f3\ = = dn. The parameters of the gauge 
transformation at second order can be obtained from the second order components of the metric 



jpi+diLJp = u^i+diU^ + ^-di [2(a[ + 2aHa\) - a\ - 4aHa>i] -^a' l diai-dia2+di(3' 2 +d' 2 i , (82) 
and 



1p ij = lij + diaxdjax - ^Sij(d k ai) 2 + 2 ydidj - ^<%d 2 ^ (h + did 2 j + djd 2 , 



(83) 



In the previous equations we used that uj^ = — 2£q /(5aH) = a\ at first order. Imposing that 7p is 
transverse (xp = Xp = in the notation of [2]) allows to solve for (3 2 and d 2 f. 



02 

d 2 i 



25a 2 H 2 
4 



dr 



2ha 2 H 2 

Imposing wp = in eq. ([82]) gives 



d- 2 [d 2 CodiCo ~ d^didj^CodjCo)] 



a\ + a 2 = uj^ + (3 2 



(84) 
(85) 

(86) 
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At second order, the other components of the metric transform as [2] 
$P = $ c + ^ fa" + 5a.ffc/ 1 + ^a 2 H 2 a^j + af + a' 2 + aiTa 2 , 



(87) 



-a 2 if 2 ai - aHa^ + 2 ($' ? + 2aH^ c ) 



2 

which, with the conditions ([511) . ([85]) and ([86]) above, give 



^(5jai) 2 - a#a 2 - ;j<9 2 /3 2 , 



$p = -|co + ^[Co+5" 2 (9iCo) 2 -3a- 4 ^(WiCo)] 

+ 175Q 6 2g2 g' 2 [2(5 l 9 J Co) 2 + 5(9 2 Co) 2 + 7^Co^5 2 Co] , (89) 



3 A 9 
* P = "5 C °-25 



Co + ld~ 2 (d t Co) 2 - 2d-%d J (d i ( d J ( ) 



+ 175Q 6 2g2 d- 2 [2(ft^Co) 2 + 5(9 2 Co) 2 + 7^Co^9 2 Co] , (90) 
24 

w£i = -2^^ 2 [5 2 Co^Co-9- 2 ^(9 2 Co^Co)] , (91) 
36/1 ji(kr)\ 2 TT 



7PJJ = "f U~ feT^T ~^fc/(W;Co) • (92) 
The transverse traceless projector above can be expanded to give 

PgT (d k CodiCo) = -d- 2 [d 2 Q 5 tJ + didjQo + 2(d 2 C cWo - a^Co^Co)] . (93) 

with 

©o = -V 2 [(a 2 Co) 2 -(^Co) 2 ] • (94) 

This result can be compared with ref. [3] by taking into account that Co = —(5/3)$ at linear 
order and on super-Hubble scales. The resulting metrics coincide up to a typo in eq. (A. 29) of 
Notice that in order to perform the calculation in Poisson gauge, one has to supplement the 
second order equations obtained in [2] by the proper second order matching with the initial condition 
provided by £ [3] . In our calculation this matching is already taken into account by the fact that we 
are always using the conserved quantity £. 

6 Discussion 

In this article we have developed a new approach to nonlinear cosmological perturbations, based on 
the equivalence between cosmological perfect fluids and derivatively coupled scalar fields. In this 
approach, perturbations can be studied using an ADM perturbative expansion of the scalar field 
action, similarly to what is done for the study of non-Gaussianities during inflation. 

Using this method, we have calculated the second order metric during matter domination in 
the comoving gauge ()25j) . where the scalar fluctuations are described by the comoving curvature 



6 To obtain the correct result the right hand side of eq. (A. 29) of [3] should be multiplied by a factor —2. 
Note that in their notation a factor of 1/2 multiplies all the second order contributions, i.e., cjj = u>u + 
while in our notation there is no factor of 1/2. 
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perturbation £. The metric, which is the main result of the paper, is given by 

500 = -i + ^T^Co) 2 , (95) 



.</..,- - ~di 2Co - d'^djCo) 2 + 3a- 4 0^ fc (^Co^Co) 



4 1 



_ H d- 2 [di( d 2 (o - d^didjidjCod 2 ^)] , (96) 

gij = a 2 ex.p[2((t)]5ij + a\j , (97) 

where 

C(t) = Co - 5^d- 2 ^(d fc Co^Co) . ( 98 ) 

Co is the perturbation generated during inflation, conserved on super-Hubble scales, and where jij 
is a traceless transverse tensor, given by 



4 

Hi = - 5 



9( i 1 



3 fer 7 5a 2 F 2 



(^Co^zCo) • (99) 



Note that at early times, when all the modes are out of the Hubble radius so that we can neglect 
gradients, the above metric is just 

ds 2 = -dt 2 + a 2 (t)e 2 <°^dx 2 . (100) 

This is the form of the metric in comoving gauge after all the modes have left the Hubble radius 
during inflation. This metric is known to remain the same during the subsequent evolution of the 
Universe until the modes re-enter the Hubble radius. This shows that our metric nicely matches the 
initial conditions set by inflation. 

As shown by eq. (|98p above, we have found that C is not conserved at second order during the 
matter dominated era. Indeed, super-Hubble modes can evolve under the presence of short sub- 
Hubble perturbations. However, we have concluded that a discussion on the physical significance 
of this variation cannot be done without the computation of the relation between C and observable 
quantities, which is beyond the scope of this article. 

Our method can be generalized to compute the metric during a cosmological era dominated by 
a fluid with constant equation of state w ^ 0, like during radiation dominance, by choosing the 
appropriate w in the linear evolution equation for C, eq. (f3Tj) . and in the third order action (|48p. 
The only difference is that in this case C is not constant inside the Hubble (see eq. (f3Tj) ). so that the 
second order solution will involve an integral over conformal time. Furthermore, one can generalize 
our method to study the evolution of perturbations in a Universe filled with several coupled fluids, by 
introducing a correspondent number of derivatively coupled scalar fields in the action, with possible 
couplings between them. As a future study, this method could be used to study the generation 
of gravitational waves from second order perturbations during radiation and matter dominated era 

[13 EH]. 
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